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Abstract 

We report a measurement of the CP-odd fraction and the time-dependent CP violation in 
B° — > D* + D*~ decays, using 657 million BB events collected at the T(4S) resonance with the 
Belle detector at the KEKB asymmetric-energy e + e~ collider. We measure a CP-odd fraction 
of R± = 0.125 ± 0.043(stat) ± 0.023(syst). From the distributions of the proper-time intervals 
between a B° — > D* + D*~ decay and the other B meson in the event, we obtain evidence of CP 
violation with measured parameters Ao*+d*- = 0.15 ± 0.13(stat) ± 0.04(syst) and So*+d*- = 
-0.96 ± 0.25(stat)+°;^(syst). 

PACS numbers: 11.30.Er, 12.15.Ff, 13.25.Hw 
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In the Standard Model (SM), the irreducible complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix gives rise to CP-violation In an T(45) event, 
the time-dependent decay rate of a neutral B meson to a CP eigenstate is given by 

-|At|/T B 

V(At) = — \l + q Ssm(Am d At) + Acos(Am d At) f, (1) 

At bo L L J ' 

where q — +1 (—1) when the other B meson in the event decays as a B° (B°), At = tcp—tta& 
is the proper-time difference between the two B decays in the event, Tgo is the neutral B 
lifetime and Arrid is the mass difference between the two B° mass eigenstates. The CP- 
violating parameters are defined as 

where A is a complex observable depending on the B° and P° decay amplitudes to the 
final state and the relation between the B meson mass eigenstates and its flavor eigen- 
states. At the quark level the B° — > D* + D*~ decay is a b — > ccd transition, where the 
tree amplitude is CKM-suppressed. The contribution of penguin diagrams in this decay 
is estimated to be at the percent level If penguin corrections are neglected, the SM 
expectations for the CP parameters are Ad*+d*- = and Sd*+d*- = — Vd*+d*- sin2fa, 
where fa = a,Yg[-V cd V* b ]/[V td V t l] and r] D *+ D *- is the CP eigenvalue of D* + D*~, which is +1 
when the decay proceeds through an S or D wave, or —1 for a P wave. A large measured 
deviation from this expectation can be a sign of new physics [3j. Recently Belle reported 
a 4.1cr CP violation effect in the B° — > D + D~ decay |4j; S was found to be consistent 
with — sin 20i whereas the measured A value indicated 3.2 a direct CP violation, which 
contradicts the SM and is not confirmed by BaBar j^]. This decay contains the same weak 
phase transition as B° — > D* + D*~ , therefore a precise measurement of the latter is vital for 
a correct interpretation. The CP-violating parameters as well as the CP-odd fraction in 
B° — > D* + D*~ decays have been measured by both Belle [(3] and BaBar Q. Here we report 
a new measurement with more than four times the statistics used in {(j). 

This analysis is based on a data sample containing 657 million BB pairs, collected with 
the Belle detector at the KEKB asymmetric-energy e + e~ collider [8| operating at the T(AS) 
resonance. The T(4S') meson is produced with a Lorentz boost 0j = 0.425 nearly along the 
z axis, defined as the direction opposite to that of the positron beam. Since the B° and P° 
are approximately at rest in the T(4S') center-of-mass (CM) frame, At can be determined 



from the displacement in z between the two decay vertices, At ~ Az/(Pjc), where c is the 
speed of light. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of Csl (Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K® mesons and to identify muons 
(KLM). A detailed description of the Belle detector can be found elsewhere 9[]. Two different 
inner detector configurations were used. A first sample of 152 x 10 6 BB pairs were recorded 
with a 2.0 cm radius beampipe and a 3-layer silicon vertex detector; for the remaining 
505 x 10 6 BB pairs, a 1.5 cm radius beampipe, a 4- layer silicon detector, and a small-cell 
inner drift chamber were used [10] . 

Charged particles are reconstructed requiring the transverse (longitudinal) distance be- 
tween the track trajectory and the interaction point (IP) to be less than 2.0 (4.0) cm. Neutral 
pions are reconstructed from pairs of photons with energies above 30 MeV and with a total 
momentum in the CM system p 77 > O.lGeV/c, which are required to have an invariant 
mass in the ranee 119 MeV/c 2 < M 77 < 146 MeV/c 2 . Neutral kaons are reconstructed via 
the decay K® — y 7r + 7r~ [ill]. The it + tt~ invariant mass is required to be within ±9MeV/c 2 
of the mass [12] and is constrained in mass and fitted to a common vertex. The n + 7T~ 
vertex is required to be displaced from the IP in the direction of the pion pair momentum. 
The neutral D mesons are reconstructed in the K~tt + , K~tt + 7t°, K~7t + 7t + tt~ , K®tt + 7t~, 
K®tt + 7t~7t° and K + K~ modes, while D + decays are reconstructed in the K~tt + tt + , K®tt + , 
K®tt + tc° and K + K"tt + modes. Unless specified otherwise charge-conjugated decays are 
implied throughout. 

Charged kaons and pions are separated using a likelihood ratio, 1Zk/-k = C(K)/ (C(K) + 
£(-7r)), constructed from ACC information, CDC dE/dx and TOF measurements. Charged 
tracks in 2- prong (3- or 4-prong) vertices are reconstructed as kaons if TZr/tt > 0.1 (0.6) and 
as pions when TZx/n < 0.9. These requirements have an efficiency of 97% (85%) for kaons 
in 2-prong (3- or 4-prong) vertices and 98% for pions, respectively, with fake rates of 18% 
(14%) for kaons and 12% for pions. The invariant mass of the D candidates must be within 
±6a (3a) of the nominal value for 2-prong (3- or 4-prong) decays, where a is the width of the 
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main component of the channel-dependent D mass resolution obtained from signal Monte 
Carlo (MC) samples and ranges from 2.6MeV/c 2 to 7.5MeV/c 2 . Candidate D* + mesons 
are reconstructed in the D°ir + and D + tt° modes. 

The pions from the D* decays are referred to as slow pions because of their low momen- 
tum. Slow charged pions are constrained to originate from the point where the D trajectory 
intersects the beam profile. The mass difference, AM = \M(D*) — M(D)\ is required 
to be within ±3(2.25) MeV/c 2 of the nominal value for the D° (D + ) channel. Finally, 
two oppositely charged D* mesons are combined to form a B° candidate. Because of the 
smaller product branching fraction and the large background contribution, we do not include 
(D + 7i°) (D~ 7i°) combinations. 

The selected D meson candidates are then subjected to mass- and vertex-constrained fits 
to improve their momentum and vertex resolution. To discriminate the signal B mesons from 
background, we use the energy difference AE = E^ M — E™ m and the beam-constrained 



mass Mb c = J (-^beam) 2 — (Pb M ) 2 , where E™ m is the beam energy in the CM system and 
E™ and p^ M are the energy and momentum of the B candidate in the CM system. After 
alle the above selection requirements are applied, there are on average 1.7 B° candidates per 
event in the large signal region. This region is defined by 5.23 GeV/c 2 < M^ c < 5.30 GeV/c 2 
and -0.14 GeV < AE < 0.14 GeV. 

We choose the B candidate with the smallest value of 



£ / AM, -AM, (PPG) \ 2 | 2 / M(A)-M(A)(PDG) \ 2 



where PDG refers to the world average measurement in 12| and % denotes the two D mesons. 



The e + e~ — >■ qq (q = u,d,s, andc) background is suppressed by requiring the ratio of the 
second- to zeroth-order Fox- Wolfram moments [l^] to be less than 0.4. 

We perform an unbinned two-dimensional maximum likelihood fit to the large signal 
region in the M^ c vs. AE plane. The probability density function (PDF) used to model 
the Mb c distribution is the sum of a signal and background component. The signal PDF is 
described with a Gaussian function while the combinatorial background is modeled with an 
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The AE signal distribution is fitted with the sum of two Gaussians 
where the width and mean of the second wide Gaussian, as well as the relative fraction 
of the two Gaussians, are fixed to the MC values. The AE background distribution is 
described with a second-order polynomial. Figure [1] shows two different projections of the 
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FIG. 1: (a) The M bc distribution for \AE\ < 0.04 GeV. (b) The AE distribution for M bc > 
5.27GeV/c 2 . The solid curve shows the result of the fit while the dotted curve is the background 
contribution. 

two-dimensional distribution and fit results. We obtain 553 ± 30 signal events in the large 
signal region. In the small signal region, denned by 5.27GeV/c 2 < M bc < 5.30GeV/c 2 and 
-0.04 GeV < AE < 0.04 GeV the signal purity is 55%. 

To obtain the CP-odd fraction we perform a time-integrated angular analysis in the 
transversity basis [3]. The differential decay rate as a function of the transversity angle is 

1 dT 



T d cos 8 



tr 



^(P + P||)sin 2 8 tT +^R ± cos 2 8 t 



(4) 



where Ro t \\ and R± are the fractions of the longitudinal, transverse parallel and transverse 
perpendicular components in the transversity basis. R$ and R» are the fractions of the 
CP-even polarization while R± is the fraction of the CP-odd one. A one- dimensional fit to 
the cos#tr distribution allows the extraction of the CP-odd fraction, where 8 tr is the polar 
angle between the momentum of the charged slow pion in its mothers D* rest frame and 
the normal to the other D* decay plain. The measured distribution of cos# t r is distorted, 
in particular due to the angular resolution of the slow pion. The shapes of the CP-odd 



and CP-even polarizations are obtained 
fraction from the previous Belle analysis 



rom a signal MC sample taking the R /(Rq + R\\) 
The background shape is obtained from the fit, 



but limited to be a symmetric polynomial, i.e. a b kg • cos 2 8 tI + 1. The signal-to-background 
ratio is determined on an event-by-event basis using the M bc — AE distribution. The fit to 
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FIG. 2: The cos#tr distribution for events in the small signal region, the points with error bars 
represent data. The solid curve is the result of the fit, the dotted curve shows the background 
contribution. The CP-even and CP-odd contributions are the dot-dashed and dashed curves, 
respectively, and are visible above the dotted background curve. 



the large signal region yields 



0.125 ±0.043 



(5) 



and abkg = —0.02 ± 0.04. The fit result is shown in Fig. [2j superimposed on the cos^tr 
distribution in the small signal region. This result is compatible with previous Belle and 
BaBar measurements jfj , Q] • 

The systematic uncertainty on R± is obtained by varying the fixed parameters within 
their errors. The signal efficiency and the Rq/{Rq + R\\) parameters give rise to systematic 
uncertainties of 0.003 and 0.009, respectively. When varying the number of signal events by 
±l<j and the signal shape in Mt, c and \AE\ such that the data points in the lower tail in 
\AE\ are well described, a systematic uncertainty of 0.003 is obtained. A fast MC is used to 
estimate any possible fit bias; we find a small shift of 0.002. Tighter vertex quality cuts lead 
to a 0.013 difference in R±. Finally, a peaking background contribution of 6.6% obtained 
from the MC is added, to which we conservatively assign a CP-odd behavior, leading to a 
0.016 change in the central value. The different contributions are summed in quadrature to 



S 



yield a systematic uncertainty of 0.023 in Rj_. 

To determine the CP-violating parameters, the signal P°-meson decay vertex is recon- 
structed by fitting the momentum vector of the D meson with the beam spot profile. No 
information on the slow pions is used. After additional requirements on the number of SVD 
hits and the vertex fit quality, we obtain 511 ± 28 events. 

The tag-side decay vertex and the flavor of the tag-side B meson are obtained inclusively 
from properties of particles that are not associated with the reconstructed B° — > D* + D*~ 
decay [17|. The PDF used to describe the At distribution is: 

V{At) = J [U g V sig (At') + (1 - / sig )P bkg (At')] ■ P rcs (Ai - At')d(At'). (6) 

(7) 

The signal fraction, / s ; g and the CP-odd probability are obtained on an event-by-event 
basis, using the previous fits to the Mt, c , AE and cos 6^ distributions. The function P sig 
is given by Eq. [1] and modified event-by-event to incorporate the effect of incorrect flavor 
assignment. A dilution factor, [1 — 2/j_(cos # t r)] for S takes into account the fraction of 
the CP-odd component. We assume 5 ev en = Sodd (Awen = A dd) and define it as S (A). 
The tagging quality is parameterized by a variable r that ranges from r = (no flavor 
discrimination) to r = 1 (unambiguous flavor assignment). The data is divided into seven r 
intervals. The wrong tag fraction w r , possible tagging performance differences between B° 
and P° decays (Aw r ), and the parameters of the resolution function P res are determined 
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id- 



using a high-statistics control sample of semileptonic and hadronic b — y c decays 
However, the width of the main Gaussian component of the resolution is determined using 
a B° — > control sample. The parameters of Pbkg(At) are obtained from a fit to 

the At distribution in sideband (M bc < 5.27GeV/c 2 ) events. 

The free parameters in the fit are Ad*+d*- and iSd*+z>-; these are determined by max- 
imizing an unbinned likelihood function for all events in the large fit region. The result 
is: 

S D * +D *- = -0.96 ±0.25, 

A D * +D *- = +0.15 ±0.13, (8) 

with a statistical correlation of 11%. The significance of CP violation using the statistical 
uncertainty only is 3.4 a. Our measurements of S and A are consistent with the SM ex- 




At (ps) 



FIG. 3: Top: At distribution of well-tagged B° — s> D* + D*~ candidates (r > 0.5) for q = +1 and 
q = — 1. The shaded area is the background contribution while the solid and dashed curves are the 
superposition of the total PDFs for well-tagged q = — 1 and q = +1 events, respectively. Bottom: 
fitted raw asymmetry of the top two distributions. 



pectation for a tree-dominated b — > ccd transition. The large direct CP violation measured 
in B° — > D + D~ [4] is thus not confirmed in this b — > ccd decay mode, in agreement with 
BaBar's result We define the raw asymmetry in each At bin as (iV + — NJ)/ (N + + A r _), 
where N + (AL) is the number of observed candidates with q — +1 (— 1). Figure [3] shows the 
At distribution and the raw asymmetry for events with a good-quality tag (r > 0.5) in the 
smalls signal region. 

The systematic uncertainties on the CP- violation parameters are summarized in Table [B 
The largest contribution comes from the R± fraction, which only affects Sd*+d*-. The 
systematic uncertainty due to the signal-to-background ratio is determined by varying the 
signal yield with ±1<t, the shape parameters such that the data points in the lower tail in 
|Ai?| are well described, and the value of R± in a correlated way, as the signal purity also 
affects the angular analysis. R± is varied by 0.003, which is the systematic error in R± 
due to the signal purity and shape. The contribution of the resolution function and the 
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TABLE I: Systematic errors on the CP-violating parameters for B — > D* + D* decays. 



Source 


A v * 




s D * 


(- d*- 


CP-odd fraction R± 


±0. 


.004 


±0. 


.109 


Signal purity and shape 


±0. 


.020 


±0. 


.030 


Standard resolution function 


±0. 


.004 


+0, 
— 0. 


000 
102 


Resolution from control sample 


±0. 


.002 


±0. 


.030 


Background shape 


±0. 


.000 


±0. 


.006 


Fit bias 


±0. 


.010 


±0. 


.031 


Am d , t b o 


±0. 


.002 


±0. 


.004 


Flavor tagging 


±0. 


Oil 


±0. 


.020 


Vertex cuts 


±0. 


.003 


±0. 


,028 


At fit range 


±0. 


.010 


±0. 


,004 


Peaking background 


±0. 


,010 


+0. 
-0. 


000 
027 


Tag-side interference 


±0. 


,034 


±0. 


,007 



Total ±0.044 +! 



background shape to the systematic error is estimated by varying each parameter by ±lcr. 
Varying the resolution parameters moves Sd*+d*- further away from zero. A fast MC is 
used to estimate the bias of the CP violating parameters for the measured values. The Am^ 
and tbo parameters are varied around their world averages [ljj]. Systematic errors due to 
uncertainties in wrong tag fractions are estimated by varying the parameters wi and Awi 
in each r region by their ±la errors. The vertex quality cut is changed to £ < 125 and the 
effect is included in the table. The At fit range is changed from At < 70 ps to At < 10 ps. 
A peaking background contribution is added with no CP violation. Finally, the tag-side 
interference uncertainty is included [ljj. The different sources are added in quadrature to 
yield ±0.04 for A and to[\l for S, reducing the significance of CP violation to 3.1 a. 

We performed various cross-checks such as a fit to the CP asymmetries of the control 
sample B° D^ + D^-, which gives A = -0.02 ± 0.03(stat) and S = -0.07 ± 0.04(stat); 
these values are consistent with no CP asymmetry. The lifetime fit to the B° — > D* + D*~ 
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sample is consistent with the world average value 12]. 

In summary, we have performed new measurements of the CP-odd fraction R± = 
0.125±0.043(stat)±0.023(syst) and CP-violation parameters A D *+ D *- = 0.15±0.13(stat)± 
0.04(syst) and S D *+ D *- = -0.96 ± 0.25(stat)+J;^(syst) for the decay B° -»■ D* + D*- using 
657 x 10 6 PP events. We obtain evidence of CP violation with 3.1a significance including 
systematic uncertainties. These measurements are consistent with and supersede our previ- 
ous results jg]] . They are also in agreement with the SM prediction for b — > c tree amplitudes 
and do not confirm the large direct CP violation seen in the B° — > D + D~ decay. 
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